Disorder-induced orbital ordering in doped manganites 
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We study the effect of quenched disorder on the ordering of orbital and magnetic degrees of 
freedom in a two-dimensional, two-band double- exchange model for e g electrons coupled to Jahn- 
Teller distortions. Using a real-space Monte Carlo method, we find that disorder can induce a 
short-range ordering of the orbital degrees of freedom near 30% hole doping. The most striking 
consequence of this short range ordering is a strong increase in the low temperature resistivity. 
The real-space approach allows to analyze the spatial patterns of the charge, orbital, and magnetic 
degrees of freedom, and the correlations among them. The magnetism is inhomogeneous on the 
nanoscale in the short-range orbitally ordered state. 
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I. INTRODUCTION 

Hole-doped perovskite manganites REi^AE^MnOa 
(RE=rare earth, AE=alkaline earth) have attracted 
great attention from the condensed matter community 
over the last decade [l|. While the initial surge of re- 
search activities on these materials was triggered by the 
discovery of the colossal magnetoresistance (CMR) ef- 
fect, a rich variety of phases and phase transitions was 
subsequently uncovered @, 3 ■ K is now widely accepted 
that the interplay among charge, spin, orbital and lattice 
degrees of freedom is the underlying cause of the com- 
plexity and richness of the physical phenomena observed 
in manganites. Recent efforts from both experiment and 
theory have highlighted the significance of quenched dis- 
order in these materials 0, Sfi 0]- Therefore, analyzing 
the effects of disorder in mang anites has become an active 
area of research @, 0, [l(| fllf . 

Disorder is generally viewed as an agent for suppress- 
ing the ordering tendencies of the microscopic degrees of 
freedom. Experiments on the half-doped (x = 0.5) man- 
ganites show that quenched disorder indeed spoils the 
long-range ordering of the charge, orbital, and spin vari- 
ables leading, in some cases, to a short-range ordering of 
these microscopic degrees of freedom 0, [(f . The oppo- 
site effect, however, is observed in manganites near 30% 
hole doping, where an ordering of the orbital degrees of 
freedom is induced by the presence of quenched disorder 

0- 

In manganites, the average ta and the variance a 2 of 
the A-site ionic radii are known to control the single- 
particle bandwidth and the magnitude of quenched dis- 
order, respectively [13 ] . Samples with constant ta and 
varying a 2 were used in the experiments of Ref. Q with 
a combination of La, Pr, Nd, Sm, and Ca, Sr, Ba at the 
A-site, while keeping x = 0.3. An increase in the low- 
temperature resistivity by four orders of magnitude was 
attributed to the onset of orbital ordering, which was 



also evidenced from the structural changes analyzed via 
powder x-ray diffraction. Magnetism is affected strongly 
with a reduction in both, the Curie temperature Tc and 
the saturation value of the magnetization. This dop- 
ing regime is also believed to be magnetically inhomoge- 
neous, as independently inferred from NMR and neutron 
scattering experiments [III, [TH . 

Disorder has been included previously in models for 
manganites to study its influence on the long-range or- 
dered phases [1, [HI, [ljj] especially near a first-order phase 
boundary or in the vicinity of phase separation [ill . [l7| . 
The idea that quenched disorder may lead to a partial 
ordering of the orbital degrees of freedom in manganite 
models has so far remained unexplored. 

In this paper, we study a two-band double-exchange 
model with quenched disorder using a real-space Monte 
Carlo method. Disorder is modelled via random on-site 
energies selected from a given distribution. We consider 
two different types of distributions, described in detail 
in the next section. Here and below we refer to these 
distributions as (i) binary disorder and (ii) random scat- 
terers. While the binary disorder has no significant effect 
on the orbital degrees of freedom, random scatterers lead 
to orbitally ordered (OO) regions, and a sharp increase 
in the low-temperature resistivity is found as observed in 
the experiment Q. The magnetic structure is inhomo- 
geneous in a restricted doping regime as observed in the 
NMR and the neutron scattering experiments. Within 
clusters, staggered orbital ordering is accompanied by fer- 
romagnetism, thus providing an example of Goodenough- 
Kanamori rules in an inhomogeneous system [l8l . [l9j . 



II. MODEL AND METHOD 

We consider a two-band model for itinerant e g elec- 
trons on a square lattice. The electrons are coupled to 
Jahn- Teller (JT) lattice distortions, f2 9 -derived 5 = 3/2 
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core spins and quenched disorder as described by the 
Hamiltonian: 

a/3 

i (ij) 

-JH^Si.ffi-A^Qi.Ti + y^Q?. (1) 

i i i 

Here, c and are annihilation and creation operators 
for e g electrons, a =|, J, is the spin index and a, f3 are 
summed over the two Mn-e s orbitals d x 2_ y 2 and d^ z 2_ r 2 1 
which are labelled (a) and (b) in what follows, t 1 ^ denote 
the hopping amplitudes between e g orbitals on nearest- 
neighbor sites and have the cubic perovskite specific 
form: C - t\ a = t, tf b = t\ h = i/3, t* b = tf a = -t/y/3, 
tab = ^ba. = t/v^i where x and y mark the spatial di- 
rections [20[ ■ The e g -electron spin is locally coupled to 
the tig spin Si via the Hund's rule coupling Jh ■ The eg- 
electron spin is given by of = YZa' ^/^Wi where 
T' 1 are the Pauli matrices. J s is the strength of the 
superexchange coupling between neighboring tig spins. 
A denotes the strength of the JT coupling between the 
distortion Qi = (Qi Xl Qi Z ) and the orbital pseudospin 

= Y^ 13 c \a<j^up c iP<j- K is & measure of the lattice 
stiffness, and we set t = 1 = K as our reference energy 
scale. 

The following two forms of on-site disorder modelling 
are used: (i) binary disorder: takes equally probable 
values ±A, (ii) random scatterers: a fraction x of the 
sites are taken to have e 4 = D, while for the other sites 
£; = 0. Although the first choice of disorder is the sim- 
plest from the model point of view, the second appears 
more realistic. In real materials a fraction x of the rare- 
earth ions is replaced by alkaline-earth ions at random 
locations. Therefore, it is likely that the disorder aris- 
ing as a consequence of this substitution is connected 
to the amount of doping. This situation is modelled by 
placing repulsive potentials on a fraction x of the sites, 
which are randomly selected. A typical measure of the 
strength of a disorder distribution is its variance. For the 
binary distribution the variance is A, while for the finite 
density x of scatterers with potential strength D, it is 
D^Jx{l — x). These two models for disorder were pre- 
viously employed in a study of half-doped manganites 
PH - The JT distortions and the t 2g derived core spins 
are treated as classical variables, and we set |Si| = 1. 
Guided by earlier estimates for the JT coupling strength 
in manganites, we fix A = 1.5 pll |. and explore the vari- 
ation in the parameters A, D, and J s . 

We further adopt the simplifying limit Jh >> t, which 
is justified and frequently used in the context of mangan- 
ites 0, [2(| ■ In this limit the electronic spin at site i is 
tied to the orientation of the core spin S;. Transforming 
the fermionic operators to this local spin reference frame 
leads to the following effectively 'spinless' model for the 



e g electrons: 

(ij) « (ij) 

-A^Q^ + y^Q?. (2) 

i i 

The new hopping amplitudes t have an additional de- 
pendence on the core-spin configurations and are given 
by: 

h!l = C os I cos % + sin % sin % . (3) 

t a p 2 2 2 2 

Here, Qi and <f>i denote polar and azimuthal angles for the 
spin Si . From now on the operator Ci a (cj Q ) is associated 
with annihilating (creating) an electron at site i in the 
orbital a with spin parallel to Sj. 

The model given by Eq. (2) is bilinear in the elec- 
tronic operators and does not encounter the problem of 
an exponentially growing Hilbert space, since all many- 
particle states can be constructed from Slater determi- 
nants of the single-particle states. The difficulty, how- 
ever, arises from the large phase space in the classical 
variables Q and S. Exact diagonalization based Monte 
Carlo (ED-MC) is a numerically exact method to treat 
such problems, and has been used extensively in the past 
[§, 0, HH. The classical variables are sampled by the 
Metropolis algorithm, which requires the exact eigenen- 
ergy spectrum. Therefore iterative ED of the Hamilto- 
nian is needed, which leads to N A scaling of the required 
cpu time, N is the number of lattice sites. The N 4 scal- 
ing makes this method very restrictive in terms of the 
achievable lattice sizes, with the typical size in previous 
studies being ~ 100 sites. Since a study of larger lattices 
is essential for analyzing the nature of inhomogencitics in 
manganite models, several attempts have been made to 
devise accurate approximate schemes 0, [24|, HH . In the 
present study we employ the travelling cluster approxi- 
mation (TCA) [25[, which indeed has been ver y su ccessful 
in analyzing similar models in the recent past jl0l.[Tll.[26|. 

III. RESULTS AND DISCUSSION 

We begin with results for bulk quantities describing the 
ordering of the magnetic and the lattice degrees of free- 
dom. We focus on the 30% hole-doped system (x = 0.3) 
for a close correspondence to the experiments in Ref. 0] • 
Fig. 1(a) shows the effect of binary disorder on the tem- 
perature dependence of the magnetization m, defined via 
to 2 = ((A -1 Si) 2 )<M). Here and below (...) av denotes 
the average over thermal equilibrium configutations and 
additionally over realizations of quenched disorder. Re- 
sults for disordered systems are averaged over 4 — 6 re- 
alizations of disorder. Clearly, the magnetism is not af- 
fected much by the presence of weak binary disorder. 
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FIG. 1: (Color online) Temperature dependence of (a) the 
magnetization m and (b) the lattice structure factor Dq(q) 
at q — (tt, tt) for various values of the disorder strength A. 
(c) and (d) show the same quantities as in (a) and (b), respec- 
tively, if the on-site disorder is modelled by random scatterers 
of strength D. The concentration of scatterers is equal to the 
hole density x = 0.3. Note the order of magnitude difference 
in the magnitudes for Dg(qo) between panels (b) and (d). 
All results are at A = 1.5 and J s = 0.05. 



This is in agreement with previous studies, which find 
that the reduction in Tq is proportional to A 2 for weak 
disorder [13, [H, HI] . 

Fig. 1(b) shows the temperature dependence of the 
q = qo = {it, tt) component of the lattice structure fac- 
tor D g (q) = i\r 2 E«<Q< ■ Qj)av e^*-*). £» Q (q ) is 
a measure for the staggered distortion order in the sys- 
tem. The lattice ordering leads to orbital ordering via the 
JT coupling. An increase with A in the low-temperature 
value of Dq (qo) suggests the appearance of orbital or- 
der. However, this effect is too weak to explain the ex- 
perimental resistivity data 0]. Moreover, the increase at 
low T in -Dq (qo) is not monotonic, which becomes clear 
by comparing the results for A = 0.4, 0.8 and 1.0 in Fig. 
1(b). 

Now we explore the results for the disorder arising from 
random scatterers of strength D. Since the disorder orig- 
inates from the replacement of RE 3+ by AE 2+ ions, the 
density of random scatterers is kept equal to the doping 
concentration x. m(T), shown in Fig. 1(c), is affected 
strongly upon increasing D, with a decrease in the sat- 
uration value of the magnetization pointing towards a 
magnetically inhomogeneous groundstate. More impor- 
tantly, a monotonic increase with D is observed in the 
low-tempcrature values of Dg(qo) (see Fig. 1(d)). The 
rise in Dq (qo) clearly indicates the emergence of orbital 
ordering in the system, with the area and/or strength of 
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FIG. 2: (Color online) Temperature dependent resistivity 
p{T) (in units of h/ne 2 ) for varying strength of (a) binary 
disorder A, and (b) random scatterers D. Note the logarith- 
mic scale in (b). (c)-(d) Low-temperature density of states 
for the two types of disorder. 



the ordered regions increasing with increasing D. 

It is expected that these orbital ordering correlations 
are reflected in the transport properties. We therefore 
compute the dc resistivity p approximated by the in- 
verse of low-frequency optical conductivity, which is cal- 
culated by using exact eigenstates and -energies in the 
Kubo-Greenwood formula [3fJ. Fig. 2(a)-(b) show p as a 
function of temperature for the two disorder models de- 
scribed above. The low-temperature resistivity increases 
upon increasing the binary disorder strength A (see Fig. 
1(a)). For small values of A the resistivity curves ap- 
pear parallel to each other below T ~ 0.1. The resistiv- 
ity therefore follows Mathiessen's rule, i.e. p{T) for the 
disordered system is obtained from p(T) for the clean 
system by simply adding a constant contribution arising 
from the scattering off the disorder potential, dp/ dT re- 
mains positive at low temperature indicating a metallic 
behavior. This oversimplified description, however, does 
not take into account the disorder induced changes in the 
orbital ordering correlations and the related changes in 
the density of states (discussed below). 

Random scatterers lead to a drastically different be- 
havior. The low-temperature rise in p(T) covers several 
orders of magnitude (see Fig. 2(b)). The negative sign 
of dp/dT for D > 1 signals an insulating behavior. Upon 
increasing the disorder strength D we therefore observe 
a metal to insulator transition. For x = 0.3 both disor- 
der models have the same variance, if A ~ 0.461? holds. 
Comparing, therefore, the results for D = 2 and A = 1, 
we have to conclude that the drastic rise in the resis- 
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FIG. 3: (Color online) Low-temperature distribution func- 
tions generated from the Monte Carlo data for the nearest- 
neighbor correlations Cq of the lattice distortions for (a) bi- 
nary disorder and (b) for random scatterers. The curves for 
different A are off-set along the y-axis for clarity. Cq is posi- 
tive (negative) for ferro- (antiferro)-distortive patterns of the 
lattice variables. 



tivity for random scatterers can not be attributed to the 
strength of the disorder potential. In fact, a large increase 
in the low-T resistivity was one of the experimental indi- 
cations for the onset of disorder-induced orbital ordering 
0- 

Figs. 2(c)-(d) highlight the difference between the 
densities of states (DOS) for the two choices of dis- 
order modelling. The DOS is defined as N(ui) = 
(N^ 1 J2i — Ei)) av , where £!j denotes the eigenvalues 
of the Hamiltonian. Here we approximate the ^-function 
by a Lorentzian with width 7 = 0.04: 



y/tt 



[Y 2 + (w - Bk 



(4) 



The DOS for the clean system has a pseudogap struc- 
ture near the chemical potential. For binary disorder, 
the pseudogap slowly fills up with increasing A. In con- 
trast, it deepens upon adding random scattering centers 
and even leads to a clean gap for D > 3. This opposite 
behavior is partly responsible for the drastically different 
low-tempcraturc resistivity discussed above. The three- 
peak structure for large values of D in Fig. 2(d) can be 
understood as follows: A fraction 2x of electronic states 
split off and form a narrow impurity band centered at an 
energy D above the Fermi level of the undoped system. 
The lower band now contains a fraction 2(1 — x) of the 
states with the Fermi level located in the middle of the 
band. This leads to a situation similar to the undoped 
system, and an energy gap originating from staggered 
orbital ordering opens at the Fermi level. 

To gain further insight into the nature of the states 
in the presence of the two types of disorder, we plot the 
distribution functions for the lattice variables in Fig. 3. 
Panel (a) shows the distribution of the nearest neighbor 
lattice correlations, Cq(i) = (l/4)X}<$Qi ' Qi+s for bi- 
nary disorder; here S denotes the four nearest neighbor 




FIG. 4: (Color online) Real-space patterns of the disorder 
potential e(i), charge density n(i), and lattice correlations 
Cq(i). Top row: binary disorder with A = 1, bottom row: 
random scatterers with D = 2. The patterns in both cases 
are shown on a 24 x 24 lattice for a single disorder realization 
at T = 0.01 and x = 0.3. 



sites of site i. A negative value of Cq(i) indicates an an- 
tiferro pattern of JT distortions, and hence a pattern of 
staggered orbital ordering. The distribution function for 
C Q is defined as P(C Q ) = {N' 1 ^ 6(C Q -C Q (i))) av ; the 
(^-function is again approximated by a Lorentzian with 
width ~ 0.04. A peak in P{Cq) centered near Cq = 0.8 
for A = indicates that the clean system has weak ferro- 
distortive/ferro-orbital correlations. Tails going down to 
Cq ~ —1.4 arise in the distribution function upon in- 
cluding binary disorder. 

The distribution function P(Cq) for random scatter- 
ers looks qualitatively different. We recall that the 
strengths of the two types of disorder are related via 
A ~ 0.46D. The low-tcmpcraturc distributions P(Cq) 
are plotted in Fig. 3(b) for random scatterers. A qual- 
itative change in the shape of the distribution function 
occurs for D = 2, where a second peak centered around 
Cq ~ — 1.2 emerges. This is a direct indication that a sig- 
nificant fraction of the system becomes orbitally ordered. 
This correlates perfectly with the strong rise in _Dg(q ) 
at low temperatures (see Fig. 1(b)), and the anomalous 
increase in the resistivity (see Fig. 2(b)). 

A real-space picture for the emergence of orbital or- 
dering is presented in Fig. 4, which displays the disorder 
potential q, the electronic density rij, and the lattice cor- 
relations Cq(i). The top row for binary disorder shows 
that the charge density closely follows the disorder po- 
tential. The local lattice correlations are centered around 
Cq = 0, which is also evident from the peak in the dis- 
tribution P(Cq) shown in Fig. 3(a). The bottom row 
in Fig. 4 shows the corresponding results for the disor- 
der potential arising from random scatterers. Since the 
doping concentration in this case coincides with the con- 
centration of the scatterers, the holes are trapped at the 
impurity sites. This leaves the surrounding effectively 
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FIG. 5: (Color online) Doping evolution of the local charge 
density m and the local lattice correlations C'q (i) for a single 
realization of random scatterers with strength D = 1, on a 
24 x 24 lattice. 



undoped and thereby induces orbital ordering. This is 
apparent from the spread of the dark-blue regions and 
their cross correlation with the charge density distribu- 
tion in the bottom row of Fig. 4. Such a picture with 
orbitally ordered regions coexisting with orbitally disor- 
dered patches describes perfectly the double peak struc- 
ture of the distribution function in Fig. 3(b). 

Although we are primarily interested in the experimen- 
tally relevant case x = 0.3, it is useful to see how the 
real-space patterns evolve as one moves from low to high 
hole densities. The undoped system is an orbitally or- 
dered insulator, which turns into an orbitally disordered 
metal upon doping [3j|. We show real-space patterns 
at three different doping concentrations in Fig. 5. The 
density of random scatterers is kept equal to the dop- 
ing fraction x. The charge density distribution is largely 
controlled by the disorder distribution. At low doping, 
disconnected orbitally disordered regions are essentially 
tied to the trapped holes. With increasing x the orbitally 
disordered regions begin to connect in one-dimensional 
snake-like patterns. By further increasing the doping and 
the concentration of scattering centers the orbitally disor- 
dered regions grow. The phenomenon of disorder-induced 
orbital ordering is likely to be present only in a narrow 
doping range near and above 0.25, because for x < 0.25, 
the system is orbitally ordered even in the clean limit. 

It is worthwhile to point to a similarity between the 
effects of disorder in the present study and in a model 
analysis for d-wave superconductors with non-magnetic 
impurities. In Ref. [32| it was found that the impuri- 
ties nucleate antiferromagnetism in their near vicinity. 
Upon increasing the impurity concentration static anti- 
ferromagnetism is observed. There seems to be a perfect 
analogy between the two situations, if one interchanges 
antiferromagnetism by orbital-ordering; both are order- 
ing phenomena with the staggered ordering wavevector 




(b) 



- 0.3 



0.2 



D Q (7t,7t) 
0.1 



0.05 



0.1 



FIG. 6: (Color online) Temperature dependence of the magne- 
tization m and the staggered lattice structure factor Dq(tt, n) 
for varying superexchange coupling strength J s . The results 
are for random scatterers with strength D = 4. 



(■7T, 7r). The (tv, tt) ordering phenomena is partially trig- 
gered by the charge inhomogencitics in both cases. An 
additional complication in the present case arises from 
the spin degrees of freedom in addition to the orbital 
variables and from the anisotropy of the hopping param- 
eters. 

As inferred above from the results for the temperature 
dependent magnetization m(T), the magnetic ground- 
state appears to be homogeneous for binary disorder, but 
may be inhomogencous in the case of doped scatterers 
(see Figs. l(a),(c)). Since the magnetism is partially 
controlled by the antifcrromagnctic superexchange cou- 
pling J Sl we study the effect of increasing J s for a fixed 
large disorder strength of random scatterers. Fig. 6(a) 
shows the result for m(T) and Fig. 6(b) the result for 
the temperature dependence of Dg(qo). The saturation 
value of m(T) as well as the onset scale for ferromag- 
netism decrease with increasing J s . More importantly, 
Z?g(q ) at low temperatures increases with increasing J s 
indicating an enhancement in the orbital ordering. For 
a homogeneous system this would mean that orbital or- 
dering and antiferromagnetism are both enhanced with 
increasing J s . This is a contradiction to the Goodenough- 
Kanamori rules, which state that an orbitally antiferro 
system should be magnetically ferro. The contradiction 
is resolved by analyzing the microscopic details of this 
complicated state providing an example where the real- 
space structures are essential for a comprehensive under- 
standing. 

We show in Fig. 7, the effect of the superexchange cou- 
pling on the real-space patterns of lattice and spin vari- 
ables. The lattice correlations are shown in the top row 
and the analogously defined spin correlations Cs(i) = 
(1/4) Si ■ S i+S in the bottom row. For J s = 0.02 
the system contains orbitally ordered nanoscale regions, 
but magnetically it appears homogeneous. For J s = 0.06 
the area of the orbitally ordered regions is enlarged and 
magnetic inhomogencitics appear. The orbitally ordered 
regions remain ferromagnetic, while the orbitally disor- 
dered regions become antiferromagnetic upon increasing 
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FIG. 7: (Color online) Real-space patterns for the lattice 
correlations Cq(i), and the analogously defined spin correla- 
tions Cs(i) for varying superexchange coupling J a . The pat- 
terns are shown for a single disorder realization at T = 0.01. 
Orbitally ordered regions tend to maintain ferromagnetism, 
while the orbitally disordered regions are more susceptible 
towards antiferromagnetism with increasing J a . 

J s . In the orbitally ordered clusters of this inhomo- 
geneous system in the selected parameters regime the 
Goodcnough-Kanamori rules are therefore fulfilled. How- 



ever, upon increasing J s further to 0.1, the antifcrromag- 
netic regions start to extend also into the orbitally or- 
dered clusters. The charge density patterns (not shown 
here) are insensitive to the increase in J s . 



IV. CONCLUSIONS 

Our analysis for a two-band double-exchange model 
for manganitcs leads us to conclude that the disorder 
induced orbital ordering in manganites near x = 0.3 is 
properly described, if the density of scattering centers 
tracks the hole concentration. Within this specific model 
of quenched disorder the induced staggered orbital order- 
ing is responsible for the orders of magnitude increase in 
the low-temperature resistivity as observed in the exper- 
iments in Ref. Q- 
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